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Abstract. A semiclassical Monte Carlo approach is adopted to study the multivalley 
spin depolarization of drifting electrons in a doped n-type GaAs bulk semiconductor, 
in a wide range of lattice temperature (40 < Tj, < 300 K) and doping density 
(10 13 < n < 10 16 cm~ 3 ). The decay of the initial non-equilibrium spin polarization 
of the conduction electrons is investigated as a function of the amplitude of the driving 
static electric held, ranging between 0.1 and 6 kV/cm, by considering the spin dynamics 
of electrons in both the T and the upper valleys of the semiconductor. Doping density 
considerably affects spin relaxation at low temperature and weak intensity of the 
driving electric held. At high values of the electric field, the strong spin-orbit coupling 
of electrons in the L-valleys significantly reduces the average spin polarization lifetime, 
but, unexpectedly, for field amplitudes greater than 2.5 kV/cm, the spin lifetime 
increases with the lattice temperature. Our numerical hndings are validated by a 
good agreement with the available experimental results and with calculations recently 
obtained by a different theoretical approach. 
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1. Introduction 

The processing of a high volume of information and world wide communication is, at 
the present, based on semiconductor technology, whereas information storage devices 
rely on multilayers of magnetic metals and insulators. Semiconductor spintronics offers 
a possible direction of technological research towards the development of hybrid devices 
that could perform logic operations, communication and storage, within the same 
material technology: information could be stored in a system of polarized electron 
spins [H|2l|3llll[5l|6l[7l[8l[9l [101 ITT] , transferred as attached to mobile carriers and 
finally detected. The possibility of obtaining long spin relaxation times or spin diffusion 
lengths in electronic materials makes spintronics a viable prospective technology. 
Nevertheless, the designers of spin devices have to worry about the loss of spin 
polarization (spin coherence) before, during and after the necessary manipulations. In 
particular, efficient injection, transport, control and detection of spin polarization must 
be carefully treated [2j. Electron-spin states depolarize by scattering with imperfections 
or elementary excitations of the medium, such as phonons. Furthermore, miniaturization 
process brings the system to experience very intense electric fields, even when the applied 
voltages are very low. This means that, for the operability of prospective spintronic 
devices, the features of spin relaxation at relatively high electric fields should be firstly 
understood. In recent years there was a proliferation of experimental works in which 
the influence of transport conditions on relaxation of spins in semiconductors has been 
investigated p2J [131 ttH EEl EEl [El fT8] . All these works are focused on the study 
of coherent spin transport at low temperatures (Tl < 30 K) and under the influence 
of weak electric fields (F < 0.1 kV/cm), except for few works [121 EE] m which spin 
depolarization has been investigated with driving fields up to 6 kV/cm. 

The temporal evolution of the spin and the evolution of the momentum of an 
electron cannot be separated. The spin depolarization rates are functionals of the 
electron distribution function in momentum space which continuously evolves with time 
when an electric field is applied to drive the transport. Thus, the dephasing rate is a 
dynamic variable that needs to be treated self-consistently in step with the dynamic 
evolution of the electron's momentum. A way to solve this problem is to describe the 
transport of spin polarization by making use of Boltzmann-like kinetic equations. This 
can be done within the density matrix approach [19] , methods of nonequilibrium Green's 
functions, as the microscopic kinetic spin Bloch equation approach [20 1 I2T ] 122 1 123 1 l2l | ITT] . 
or Wigner functions [25] [26] , where spin property is accounted for starting from quantum 
mechanics equations. 

Another way is to use a semiclassical Monte Carlo approach, by taking into 
account the spin polarization dynamics with the inclusion in the code of the precession 
mechanism of the spin polarization vector [27] [28] EHJ [30] EU [32] [331 El]. 

Both methods allow to include the relevant spin relaxation phenomena for electron 
systems and take into account the details of electron scattering mechanisms, material 
properties and specific device design; their predictions have been demonstrated to be in 
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good agreement with experiments. 

Theoretical descriptions of the transport of spin-polarized electrons have been also 
achieved by the drift-diffusion approximation. The existing drift-diffusion schemes can 
be classified into two approaches accounting the spin degree of freedom differently: 
the two-component drift-diffusion model and the density-matrix based approximations 
[351 |36| [37] . General conditions for the applicability of these approximations are not 
different from the usual conditions of applicability of drift-diffusion approximations. 

Despite decades of studies, most of theoretical or simulative works have considered 
only the central valley T since the spin-orbit coupling parameters of the upper conduction 
bands have been only recently theoretically calculated by Fu et al. [HE]- Monte Carlo 
approaches have been widely adopted by groups of scientists to study spin polarized 
transport in 2D channels, heterostructures, quantum wells, quantum wires [271 EEl EQl 
I3T1 132| 133] . However, till today, to the best of our knowledge, in semiconductor bulk 
structures a theoretical investigation of the influence of transport conditions on the spin 
depolarization in the presence of high electric fields, comprehensive of the effects of both 
lattice temperature and impurity density, is still lacking. 

Aim of this work is to numerically estimate the spin lifetimes of an ensemble 
of initially-polarized electrons drifting in doped n-type GaAs bulks, with the lattice 
temperature, doping density and electric field amplitude ranging in a wide interval of 
values, focusing on the effects due to the inclusion of upper valleys of the semiconductor 
on the mechanism of depolarization. 

The paper is organized as follows: in Sec. [2] we briefly describe the multivalley model, 
used for the study of the spin depolarization dynamics and the Monte Carlo simulator; 
in Sec. |3]the numerical results are given and discussed. Final comments and conclusions 
are given in Sec. HI 

2. Theory and Monte Carlo approach 

2.1. Spin dynamics and multivalley model 

Spin dephasing may be caused by interactions with local magnetic fields originating 
from nuclei and spin-orbit interactions or magnetic impurities. The most relevant spin 
relaxation mechanisms for an electron system under non degenerate regime are: (i) 
the Elliott- Yafet (EY) mechanism, in which electron spins have a small chance to flip 
during each scattering, due to the spin mixing in the conduction band [391 00]; (h) the 
Dyakonov-Perel (DP) mechanism, based on the spin-orbit splitting of the conduction 
band in non-centrosymmetric semiconductors, in which the electron spins decay due to 
their precession around the k-dependent spin-orbit fields (inhomogeneous broadening) 
during the free flight between two successive scattering events [^Tl W2\ [7J ; (iii) the Bir- 
Aronov-Pikus (BAP) mechanism, in which electrons exchange their spins with holes 
[13]. Hyperfine interaction is another mechanism, usually important for spin relaxation 
of localized electrons, but ineffective in metallic regime where most of the carriers are 
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in extended states 

Previous theoretical [2H [XT] and experimental [UJ investigations indicate that the 
the EY mechanism is totally irrelevant on electron spin relaxation in n-type III-V 
semiconductors. Hence, in this work we analyze the spin depolarization of drifting 
electrons in n-type GaAs semiconductors by considering only the D'yakonov-Perel 
process. 

By following the semiclassical formalism, the term of the single electron Hamiltonian 
which accounts for the spin-orbit interaction can be written as 

H so = ^a-n. (1) 

It represents the energy of electron spins precessing around an effective magnetic field 
[B = hQ/fisg] with angular frequency Q, which depends on the orientation of the 
electron momentum vector with respect to the crystal axes. Near the bottom of each 
valley, the precession vector can be written as |SJ EE] 

fir = MUK ~ k 2 z )x + k y (k 2 z - kl)y + k z (k 2 x - k 2 y )z] (2) 

in the T- valley, and 

tt L = -j^x(k y - k z ) + y(k z - k x ) + z(k x - k y )] (3) 

in the L- valleys, located along the [111] direction of the crystallographic axes. In 
equations ()2])-(j3]), ki {i = x, y, z) are the components of the electron wave vector. (3r and 
(3l are the spin-orbit coupling coefficients, crucial parameters for the simulation of spin 
polarization. Here, we assume /3^=0.26 eV/A-2/H, as recently theoretically estimated 
[38]. In T- valley we consider the effects of nonparabolicity on the spin-orbit splitting by 
using [35] , 

= ah 2 ( _ E(k)9-7 V + 2rj i \ 

m^hn~E g \ E g 3-7? J 1 ' 

where a = 0.029 is a dimensionless material-specific parameter, rj = A/(E g + A), with 
A = 0.341 eV the spin-orbit splitting of the valence band, E g is the energy separation 
between the conduction band and valence band at the F point, m the effective mass and 
E(k) the electron energy. 

The quantum-mechanical description of electron spin evolution is equivalent to that 
of the classical momentum S experiencing the effective magnetic field, as described by 
the equation of motion 

^ = n x s. (5) 

dt v ; 

Every scattering event changes the orientation of the effective magnetic field B (that 
strongly depends on k) and the direction of the spin precession axis. 
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2.2. The Monte Carlo algorithm 

The Monte Carlo approach is one of the most powerful methods to simulate the transport 
properties in semiconductor devices beyond the quasi-equilibrium approximations 
[HJ HU [50] . In fact, owing to its flexibility, the Monte Carlo method presents the 
remarkable advantage of giving a detailed description of the particle motion in the 
semiconductor by taking into account the details of collisions with impurities, phonons, 
etc., specific device design and material parameters, and allows us to obtain all the 
needed information, such as the average velocity of electrons, temperature, current 
density, etc., directly without the need of calculating first the electron distribution 
function. In simulations, between two successive scattering events, each electron 
propagates along a classical trajectory and, according to the classical equations of 
motion, it is affected by the presence of external fields. The time interval between two 
collisions (time of free flight), the scattering mechanisms, the collisional angle, and all 
the parameters of the problem are chosen in a stochastic way, making a mapping between 
the probability density of the given microscopic process and a uniform distribution of 
random numbers. 

In our code the conduction bands of GaAs are represented by the T-valley and 
by four equivalent L valleys. We do not consider the X-valleys because, even for the 
highest value of the driving field (F= 6 kV/cm), the percentage of electrons in these 
valleys is always lower than 0.1%. The algorithm includes: (i) the intravalley scattering 
with acoustic phonons, ionized impurities, acoustic piezoelectric phonons, polar optical 
phonons, and for the the L-valleys also the scattering with optical nonpolar phonons; 
(ii) the intervalley scattering with the optical nonpolar phonons among the two valleys. 
The complete set of n-type GaAs parameters used in our calculations is listed in Ref. 
[5Tj . The scattering probabilities are calculated by the Fermi Golden Rule and assumed 
to be both field and spin independent; accordingly, the influence of the external fields is 
only indirect through the field-modified electron velocities. Nonlinear interactions of the 
field with the lattice and bound carriers are neglected. We also neglect electron-electron 
interactions and consider electrons to be independent [27]. The spin polarization vector 
is included into the Monte Carlo algorithm and calculated for each free carrier. From 
Eq. (|5]), the Monte Carlo (MC) simulator calculates the electron spin precession, by 
taking into account the scattering-induced deviations of precession vector suffered after 
each collision. 

All simulations are performed by using a temporal step of 10 fs and an ensemble 
of 5 • 10 4 electrons to collect spin statistics. We assume that all donors are ionized and 
that the free electron concentration is equal to the doping concentration n. 
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3. Numerical Results 

3.1. Spin lifetime calculation 

The initial non-equilibrium spin polarization decays with time as the electrons, driven 
by a static electric field, move through the medium, experiencing elastic and anelastic 
collisions. Since scattering events randomize the direction of Q, during the motion, 
the polarization vector of the electron spin experiences a slow angular diffusion. The 
dephasing of each individual electron spin produces a distribution of spin states that 
results in an effective depolarization, which is calculated by ensemble-averaging over the 
spin of all the electrons. 

The simulation of the spin relaxation starts with all the electrons of the ensemble 
initially polarized ({S) = 1) along the x-axis at the injection plane (xq = 0). After 
a transient time of typically 10 4 time steps, long enough to achieve the steady-state 
transport regime, the electron spins are initialized, the spin relaxation begins and the 
quantity (S) is calculated as a function of time. In order to extract the characteristic 
time t of the spin relaxation, the obtained trend of the spin dephasing is fitted by the 
following exponentially time decaying law 

(S x )(t)=A-exp(-t/r), (6) 

with A a normalization factor. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Time (ps) 



Figure 1. Average electron spin polarization (S x ) as a function of time, by only 
considering the electrons drifting into the T- valley (dashed line), into the L-valleys 
(dotted line) and into the (r + L)-valleys (solid line). F — 5 kV/cm, n — 10 13 cm" 3 
and T L = 300 K. 

In Fig. [T]we show the average electron spin polarization (S x ), in the presence of 
a driving electric field, with amplitude F = 5 kV/cm and directed along the a;- axis, 
with density n = 10 13 cm -3 and lattice temperature Tl = 300 K. This value of field 
amplitude is high enough to allow almost 21% of all electrons to visit the L valleys. The 




Figure 2. Spin lifetime r as a function of the electric field amplitude F, at Tl=77 
K (blue lines) and K (red lines), for three values of doping density, namely 

n = 10 13 cm -3 , n = 10 15 cm~ 3 and n — 10 16 cm -3 , with the electrons drifting in both 
the r and the L- valleys. 




Figure 3. (a) Momentum scattering time t p and (b) spin precession frequency fl as 
a function of the electric field amplitude F, at Tl = 77 K (blue lines) and Tj, = 300 K 
(red lines), for three values of doping density, namely n = 10 13 , 10 15 and 10 16 cm -3 . 



curves represent the decreasing trend of (S x ) vs. time by firstly considering only the 
electrons drifting into the T- valley (dashed line), secondly, by solely taking into account 
the electrons moving into the L-valleys (dotted line) and, finally, by considering the 
electrons moving into both the F and the L- valleys (solid line). We find a significant 
reduction of the average spin polarization lifetime caused by the spin-orbit coupling in 
L-valleys stronger with respect to that in T-valley, according to the theoretical results 
obtained by Zhang et al. (23] on quantum wells. The transition of about 21% of electrons 
to the L-valleys leads to an increase of efficacy of the dephasing mechanism, which brings 
to a reduction of (S x ) over time in the range 15 -j- 20% 
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Figure 4. Spin lifetime r as a function of the electric field amplitude F, at n = 10 
cm~ 3 for different values of lattice temperature, namely T^=77, 150, 220, 300 K, with 
electrons drifting in both the T and the L-valleys. 



3. 2. Effects of temperature and doping density on spin relaxation 

In Fig. [21 we show the spin depolarization time r as a function of the electric field 
amplitude F, for two values of the lattice temperature, namely Tl = 77 (blue curves) 
and 300 K (red curves) and three values of doping density n = 10 13 , 10 15 and 10 16 
cm -3 , respectively, leaving the electrons free to drift in both the T and the L- valleys. 
Except for the case at the doping density n = 10 13 cm -3 and Tl = 77 K, we find that r is 
always a monotonic decreasing function of F. In fact, when the field amplitude becomes 
larger, the electron momentum k increases, causing a stronger effective magnetic field, as 
expected by Equations (J2])-(J3]). Consequently, the electron precession frequency becomes 
higher, inducing a faster spin relaxation |JT] • 

For field amplitudes greater than 2.5 kV/cm, at Tl = 300 K we find depolarization 
times longer than those obtained at Tl = 77 K. In order to avoid that the observed 
behaviour could be ascribed only to stochastic fluctuations of MC computations, we 
have calculated the statistical error associated to our simulated data. We have repeated 
our simulations ten times, finding a maximum spread of 0.05 ps, which corresponds to 
about 1% of the observed variation of r with the temperature. 

We have investigated the counterintuitive behavior of longer average spin lifetimes 
obtained for hotter electrons, by adopting the proportionality law of Perel |41j : 



where Q is the spin precession frequency and r v the momentum characteristic relaxation 
time, corresponding to the time scale of the scattering events. We have calculated the 
spin precession frequency and the momentum characteristic scattering time for each 




T 



(7) 
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electron of our ensemble and in Fig. [31 we show the average momentum scattering time 
r p (panel (a)) and the average spin precession frequency Q (panel (b)) as a function of 
the electric field amplitude F, at Tl = 77 and 300 K, for the three values of doping 
density n = 10 13 , 10 15 and 10 16 cm -3 . 

The panel (a) of Fig. [3] shows that r p is a monotonically decreasing function of F 
for every values of n and Tl- At the higher temperature (Tl = 300 K) we find low values 
of t p , because electrons experience a greater number of scattering events, both in the T- 
valley and in L-valleys. Moreover, the curves at room temperature are characterized by 
only a slight slope, because in this case the thermal energy of the electrons is dominant 
with respect to the drift kinetic energy. At Tl — 77 K and for very low values of the 
electric field amplitude, since the scattering events are mainly due to ionized impurities, 
t p is greatly dependent on n, increasing its value at lower densities. At Tl = 300 K, t p 
is nearly independent on the doping density since the dominant scattering mechanism 
is due to the optical phonons. 

The panel (b) of Fig. |3] shows that the spin precession frequency Q is an increasing 
monotonic function of F. For F < 3 kV/cm, independently from the values of the doping 
density, the values of Q obtained at room temperature are larger than those obtained 
at Tl = 77 K. The increase of the spin precession frequency for electrons moving at 
higher temperatures is explained by the increasing number of electron transitions from 
the r to L valleys, being the value of the spin-orbit coupling coefficient in the L-valleys 
one order of magnitude greater than that of T-valley At Tl = 77 K, for F < 3 
kV/cm, the percentage of electrons in the central valley V is practically unitary and 
the spin precession frequency increases as the third power of electron momentum k, 
which increases with F according to the Eq. fl2]). When F is greater than 3 kV/cm, 
the percentage of electrons in the L-valley is high enough to lead Q for having a nearly 
linear trend (see Eq. (j3J)). At Tl — 300 K, for F < 3 kV/cm, the term of thermal energy 
is dominant with respect to the drift kinetic energy and Q vs. F shows a more slight 
increase. 

For F > 3 kV/cm, independently from the values of Tl and n, the action of F wins 
on the disorder due to the lattice temperature. In fact, except for the data obtained at 
Tl = 300 K and n = 10 16 cm -3 , which show lower values of Q, all curves coincide. The 
detached curve Tl = 300 K and n = 10 16 cm -3 is ascribed to a strong reduction of the 
percentage of electrons present in the L-valleys. 

The calculation of the square of the spin precession frequency times the momentum 
relaxation time as a function of the electron energy for each electron of the ensemble 
shows that for F > 2.5 kV/cm the average value of Q 2 t p obtained at Tl = 77 K is 
greater than that at Tl = 300 K. This finding explains the longer lifetimes observed at 
higher temperatures for field amplitudes greater than 2.5 kV/cm . 

To highlight on the nonmonotonic electric field dependence of r, observed at the 
doping density n = 10 13 cm -3 , we have investigated the spin lifetime dependence on F 
also for different values of the lattice temperature, namely 7^=77, 150, 220, 300 K (see 
Fig. HJ). Up to T L = 150 K the decoherence times slightly depend on the temperature 
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and are characterized by a marked maximum. The presence of a maximum in the spin 
depolarization time can be explained by the interplay between two competing factors, 
both due to the increase of the electric field. In the momentum space, at greater field 
amplitudes, the electrons occupy states with larger k, characterized by a stronger spin- 
orbit coupling, causing an enhancement of the spin inhomogeneous broadening. On 
the other hand, a larger electric field also brings about an increase of the number of 
scattering events, giving rise to a reduction of the momentum relaxation time. This in 
turn causes an increase of the spin relaxation time as follows from Eq. UJ 

At low values of temperature and for electric field amplitudes 0.1 < F < 0.5 kV/cm 
the inhomogeneous broadening is still marginal and the spin relaxation phenomenon is 
dominated by the momentum scattering. In particular, the number of electron scattering 
events, which are mainly due to interactions with acoustic phonons at very weak electric 
fields, increases its value because of the triggering of the scattering mechanism by ionized 
impurities, causing a reduction of r p . For field amplitudes greater than ps 0.5 kV/cm, 
the enhancement of the spin-orbit coupling, which is k-cubic dependent, is faster than 
the decrease on r p . Consequently, the spin lifetime starts to decrease with the increasing 
of the electric field. The nonmonotonic electric field dependence of r is not observed 
for Tl > 150 K where, because of the greater drift electron velocities, the loss of spin 
polarization is mainly due to the strong effective magnetic field. 

3.3. Comparison with experiments and with other theoretical approaches 

Unfortunately, till today experimental investigations of the ultrafast relaxation of 
electron spin during drift transport in bulk semiconductors, at both sample temperatures 
higher than 30 K and applied field amplitudes greater than 0.1 kV/cm, are still missing. 

Although the main aim of this work is the investigation of the influence of transport 
conditions on the spin relaxation of electrons driven by high-intensity electric fields, 
including the effects of choosing both the lattice temperature and impurity density 
in a wide range of values, in order to validate the prediction capability of our MC 
code, we have performed a comparison between our numerical spin relaxation times 
and very recent experimental results on the electron-spin-relaxation rate, obtained by 
Ref. [52]. These experiments were carried out by performing spin noise spectroscopy 
on a sample of n-type GaAs at a doping concentration of n = 2.7 ■ 10 15 cm" 3 , without 
any driving field and for lattice temperatures Tl between 4 and 80 K. In Fig. [5] we plot 
the temperature dependence of the spin-relaxation rate calculated from our code (solid 
line), together with the experimental data (circles, sample B in Ref. [52]). In order to 
best fit the experimental points with our numerical trend, we have utilized the spin-orbit 
coupling coefficient in T- valley (3r as a free parameter, obtaining the best agreement with 
/?r=19 eV-A 3 . This value is only slightly different from the value (23.8 eV-A 3 ), recently 
estimated by using the tight binding theory [38]. However this value is still within the 
reasonable range of values calculated and measured via various methods, as reported in 
Ref. [53]. For sample temperatures greater than 45 K, our numerical trend well agrees 
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with the experimental data, while at lower temperatures a considerable discrepancy is 
found, probably due to the neglecting the electron-electron scattering mechanism. In 
fact, as shown in Ref. [21] at T^=40 K, the assumption of neglecting the Coulomb term 
gives spin lifetime values smaller (hence relaxation rates higher) than those obtained 
from full calculation including all scattering mechanisms. 



800 




200 - 



40 50 60 70 80 90 100 

sample temperature T L (K) 

Figure 5. Temperature depended measurements of the spin-relaxation rate from the 
experiment (sample B) in Ref. |52] (circles) and numerical data obtained from our 
Monte Carlo code (solid curve), n = 2.7 • 10 15 cm~ 3 , /3r=19 eV-A 3 

In order to further test the effectiveness of our code we have compared our one- 
valley numerical data with the calculation of the effects of a low-amplitude electric 
field (F < 2 kV/cm) on spin relaxation in n-type III-V semiconductor bulks, lately 
obtained from the fully microscopic kinetic spin Bloch equation (KSBE) approach |24j. 
To the best of our knowledge, that paper is the only one in which the electric field 
dependence of spin lifetime has been investigated. In Fig. [6] we plot the ratio of the 
spin relaxation time under electric field to the electric-field-free one t(F)/t(F = 0) and 
the ratio between the hot-electron temperature and the lattice temperature T e /T L , as 
a function of the applied field obtained from our Monte Carlo code for a GaAs bulk 
with n = 10 16 cm -3 at T^=300 K. These results are compared with the calculations 
from the KSBEs (see Fig. 15(b) of Ref.[24j). Our finding for the spin lifetime well 
agrees with the theoretical results in all the investigated range, while the values of the 
electron temperature are systematically slightly lower than those obtained from the 
KSBE approach. Once more time, this small deviation can be due to the fact that the 
Coulomb scattering is neglected in our computation. In fact, the electron-electron (e-e) 
scattering is very important in determining the hot-electron temperature and the hot- 
electron temperature influences both the electron-longitudinal optical (e-LO) phonon 
and the electron-impurity scattering. Hence the e-e scattering effectively influences 
the spin relaxation. However, since at high temperatures, the e-e scattering is weak 
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compared to e-LO phonon scattering, it leads only to a marginal decrease in the spin 
relaxation [23]. Moreover, Romer et al. [52] have experimentally shown that, for doping 
concentrations below the metal-to-insulator transition, the electron-electron interaction 
is weaker at low impurity densities. In our finding this decrease is surely negligible for 
the two lower values of donor concentrations, namely n = 10 13 cm -3 and n = 10 15 cm -3 . 

The analysis of the influence on the spin lifetime of the intravalley and intervalley 
Coulomb scattering at high lattice temperatures is a very interesting problem which 
deserves further detailed investigations that will be the subject of a forthcoming paper. 

10 , , , , 1 2 
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0.5 1 1.5 2 
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Figure 6. Ratio of the spin relaxation time under electric field to the electric-ficld- 
free one t(F)/t(F — 0) and ratio between the hot-electron temperature and the lattice 
temperature T 6 /Tl as a function of the applied field obtained from our Monte Carlo 
code, n = 10 16 cm" 3 , T L = 300 K. 



4. Conclusions 

A full understanding of the role played by the lattice temperature, the doping density 
and the amplitude of high-intensity electric field on the electron spin dynamics in 
semiconductors is essential for the design and the fabrication of spintronic devices. 

In this work we have studied the spin lifetimes of an ensemble of conduction 
electrons, drifting in doped n-type GaAs bulk crystal, by using a semiclassical Monte 
Carlo transport model. The spin depolarization is investigated in a wide range of lattice 
temperatures and doping densities by simulating the electrons driven by an electric field 
with amplitude 0.1 < F < 6 kV/cm, including both the T- valley and the L- valleys into 
the spin dephasing dynamics. 

Our results show that the electron spin lifetime is not marginally influenced by 
the driving electric field, the lattice temperature and the impurity density, which hence 
represent key parameters into the depolarization process. We find a significant reduction 
of the average spin polarization lifetime at high values of the electric field, caused by the 
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stronger spin-orbit coupling of electrons in the L-valleys. In the nondegenerate regime 
the doping density considerably affects spin lifetimes at nitrogenum temperature and 
weak intensity of the driving electric field. Moreover, for field amplitudes greater than 
2.5 kV/cm, we observe spin lifetimes longer at room lattice temperatures with respect 
to those observed at Tl = 77 K. 

The prediction capability of our numerical code has been positively validated 
through a comparison with the available experimental results and with the calculations 
obtained by a different theoretical approach. 
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